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Decoherence of a localized electron spin in a solid state material (the "central spin" problem) at 
low temperature is believed to be dominated by interactions with nuclear spins in the lattice. This 
decoherence is partially suppressed through the application of a large magnetic field that splits the 
energy levels of the electron spin and prevents depolarization. However, the dephasing decoherence 
resulting from a dynamical nuclear spin bath cannot be removed in this way. Fluctuations of the 
nuclear field lead to an uncertainty of the electron's precessional frequency in a process known as 
spectral diffusion. This paper considers the effect of the electron's wavefunction shape on spectral 
diffusion and provides wavefunction dependent decoherence time formulas for a free induction decay 
as well as spin echoes and concatenated dynamical decoupling schemes for enhancing coherence. 
We also discuss a dephasing of a qubit encoded in singlet-triplet states of a double quantum dot. 
A central theoretical result of this work is the development of a continuum approximation for the 
spectral diffusion problem which we have applied to GaAs and InAs materials specifically. 
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I. INTRODUCTION 

Understanding quantum decoherence is a fundamental 
subject of interest in modern physics. In this work, we 
theoretically study the issue of quantum decoherence for 
the problem of one localized electron spin in a solid state 
nuclear spin environment, where the electron spin even- 
tually loses its quantum phase memory (i.e., dephases) 
due to its interaction with the surrounding nuclear spin 
bath. This is often called "central spin" decoherence in a 
spin bath, with the localized electron spin being the cen- 
tral spin and the surrounding nuclear spin environment 
being the spin bath. This particular problem is impor- 
tant in the context of quantum information processing 
and quantum computation using localized electron spins 
as qubits, and as such, we concentrate on a few sys- 
tems of interest in solid state quantum computation ar- 
chitectures, namely, Si:P donor electron spin qubits and 
GaAs and InAs quantum dot electron spin qubits, all of 
which have considerable recent experimentalLSiLlSi and 
theoreticalSiLM.i2iiiiiSii2.ili5, interest. The theory we 
develop is, however, applicable to the general situation of 
the quantum dephasing of a single localized electron spin 
in solids due to the environmental influence of the slowly 
fluctuating nuclear spin bath consisting of many millions 
of surrounding nuclear spins mutually flip-flopping due 
to their magnetic dipolar coupling. 

The issue of specific interest in this paper, as the title 
of this paper suggests, is how the detailed form of the con- 
finement for the localized electron in the solid (e.g., the 
exponentially confined hydrogenic confinement for the 
localized P donor electron state in Si or the Gaussian- 
type simple harmonic oscillator confinement for the lo- 
calized electron in the GaAs/InAs quantum dot) could 
have qualitative influence on its nuclear induced spin de- 
phasing. This subtle (but potentially significant) depen- 



dence of electron spin dephasing on the nature of the 
electron localization has recently been emphasized in the 
discovery^ that a particular type of dynamical decou- 
pling (DD) sequence^ can be ideal in restoring quantum 
coherence in the GaAs quantum dot system, but not par- 
ticularly effective in the Si:P system, which can be traced 
back to the Gaussian versus the exponential wavefunc- 
tion localization in the two systems, leading to the valid- 
ity or the lack thereof of a particular time perturbation 
expansion as discussed in depth in this work. Thus, a 
detailed investigation of the effect of the localized elec- 
tron wavefunction on the nuclear induced electron spin 
dephasing problem is both important and timely in view 
of the intense current activity in fault-tolerant quantum 
computation using spin qubits in semiconductors. 

It is important in the context of studying electron spin 
decoherence to distinguish among three different spin re- 
laxation or decoherence times, Tf, T 2 *, and T2, which are 
discussed in the literature. (We should mention right 
at the outset that our work is focused entirely on T2, 
sometimes also denoted Tm or spin memory time. T2 
is variously called spin decoherence time, spin dephas- 
ing time, transverse spin relaxation time, spin-spin relax- 
ation time, and spin memory time in the literature.) The 
spin relaxation time T1 , also often called the longitudinal 
spin relaxation time or the energy relaxation time, is con- 
nected with the spin flip process which, in the presence 
of an externally applied magnetic field (the case of inter- 
est to us in this work), necessarily requires phonons (and 
spin-orbit coupling) to carry away the electron spin Zee- 
man energy, which is 3 orders of magnitude larger than 
the nuclear spin Zecman energy. This Ti-relaxation pro- 
cess can be made arbitrarily long by lowering the lattice 
temperature so that phonons are simply not available to 
provide the energy conservation. At the low (~ 100 mK 
or lower) temperatures of interest to us in the quantum 



2 



computing context, the relevant T\ times are very long 
(Ti > 100 ms 3> T 2 ) and are unimportant for our con- 
sideration. The T 2 time is the relevant decoherence time 
in the presence of substantial inhomogeneous broadening 
as, for example, in ensemble measurements over many 
electron spin qubits with varying (i.e., inhomogeneous) 
nuclear spin environments. In the context of single spin 
qubits, i.e. involving a single electron spin, the T 2 * deco- 
herence sets in due to the requisite time averaging which, 
due to ergodicity, becomes equivalent to the spatial in- 
homogeneity of the varying nuclear spin environments of 
many electron spins. Thus, T 2 * is measured cither in a 
measurement over an ensemble of localized spins with 
the associated spatial averaging or in a time-averaged 
measurement for a single spin over many runs. A spin 
echo (or Hahn spin echo) measurement gets rid of the in- 
homogeneous broadening and characterizes T 2 , the pure 
dephasing time of a single spin (typically for the systems 
of our interest T 2 * < T 2 /1000 and T 2 < Ti/1000), which 
is what we theoretically study in this work. A closely re- 
lated, but by no means identical, definition of T 2 comes 
from considering the free induction decay (FID) of a sin- 
gle spin in a single-shot measurement without involving 
either spatial averaging over many spin qubits or tem- 
poral averaging over many runs. Alternatively, FID is 
observed in a homogeneous ensemble. We will call such 
a FID dephasing time Tj (< T 2 ) to distinguish it from 
the spin echo dephasing time T 2 . 

The above discussion of Ti, T 2 , T 2 , and Tj illustrates 
the considerable semantic danger of discussing "spin de- 
coherence" because, depending on the context, the "spin 
decoherence time" for the same system could vary by 
many orders of magnitude (i.e., T\ > T 2 > Tj > T 2 * , 
etc.). To avoid such confusion, we emphasize that, in 
our opinion, the only sensible way of discussing spin deco- 
herence is by considering specific experimental contexts. 
Our definition of T 2 is thus the decoherence time mea- 
sured in a Hahn spin echo experiment. The only license 
we take with our definition of T 2 is that we continue using 
T 2 as the notation for spin decoherenc time even in situa- 
tions where the spin coherence has been extended far be- 
yond the Hahn spin echo time by using multiple pulse se- 
quences [Carr-Purcell-Meiboom-Gill (CPMG), concate- 
nated dynamical decoupling (CDD), etc.] much more 
complex than the single 7r-pulse Hahn sequence. For us, 
therefore, T 2 is the spin decoherence time as measured in 
an echo-type pulse sequence measurement, which could 
be a simple 7r-pulse spin echo or more complex pulse se- 
quences meant to prolong spin memory beyond the spin 
echo refocusing. 

Finally, we point out an additional important compli- 
cation, often erroneously neglected in the literature, as- 
sociated with discussing spin decoherence in terms of a 
single decoherence time parameter, T Coh (e.g., T\ or T 2 
or T 2 or Tj, etc.). Such a description assumes, by def- 
inition, that the quantum memory (i.e., some precisely 
defined quantum amplitude or probability) falls off in a 
simple exponential manner with time, i.e., exp (— t/T Coh ) 



or exp (— [t/T Coh ] n ), where n is a constant, so that a sin- 
gle decoherence time T Cah can completely parametrize the 
nature of decoherence. This is, however, not always the 
case, and the detailed functional dependence of quan- 
tum coherence on time almost always changes with t in a 
complex manner, ruling out any simple single-parameter 
characterization of spin decoherence. To be consistent 
with the standard literature, we often discuss or describe 
our results by a single T 2 , but we simply define T 2 as 
the time it takes for the quantum memory to decay by a 
factor of e (or the extrapolated time at which an approx- 
imate exponential decay form will reach 1/e). This way 
we are not assuming any particular functional form of the 
quantum memory versus time decay. To be explicit, our 
results clearly indicate the decay of the spin probability 
density over time. 

The rest of the paper is organized as follows: SectionlTTl 
introduces the concept of spectral diffusion, which is the 
only spin dephasing mechanism considered in this work 
(we believe it to be the most important spin decoherence 
mechanism for solid state quantum information process- 
ing using electron spins). Sections III Al III Bl and [TlC] for- 
mally define the problem in terms of the Hamiltonian, the 
decoherence measure, and pulse sequences, respectively. 
In Sec. IIIII we review our cluster expansion method^ ' 11 ! 15 
for solving this problem. Section IIVI introduces the role 
of confinement or wavefunction by considering the initial 
time ( "short time" ) decay of spin coherence and then 
providing the detailed theoretical considerations associ- 
ated with the functional form of the localized wavefunc- 
tion as relevant for spin dephasing; Section IIVCI con- 
tains a particularly important continuum approximation, 
which provides convenient formulas that yield estimated 
T 2 times as a function of the wavefunction size and shape. 
In Sec. El we consider a specific recent experimental sit- 
uation of singlet-triplet states in a double quantum dot 
and show its equivalence to the single electron case as far 
as spin dephasing is concerned; Sec. IV Al discusses the 
limit on the experimentally discovered Zamboni effect in 
enhancing spin coherence. In Sec. IVIi we conclude with 
a summary and a brief discussion of open questions. 



II. SPECTRAL DIFFUSION 

The spin decoherence mechanism known as spectral 
diffusion (SD) has a long history^^^^^^, and has 
been extensively studied recently 1 ^ Wi 8 i 9 i u i 12 i 13 i 14 i 15 i 24 
in the context of spin qubit decoherence. Consider a lo- 
calized electron in a solid. This is our central spin. The 
electron spin could decohere through a number of mech- 
anisms. In particular, spin relaxation would occur via 
phonon or impurity scattering in the presence of spin- 
orbit coupling, but these relaxation processes are strongly 
suppressed in localized systems and can be arbitrarily re- 
duced by lowering the temperature and applying a strong 
external magnetic field, creating a large electronic Zee- 
man splitting. In the dilute doping regime of interest in 
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Spectral diffusion of a Si:P spin 




FIG. 1: (Color online) The electron of a P donor in Si ex- 
periences spectral diffusion due to the spin dynamics of the 
enveloped bath of Si nuclei. Of the naturally occurring iso- 
topes of Si, only 29 Si has a net nuclear spin, which may con- 
tribute to spectral diffusion by flip-flopping with nearby 29 Si. 
Natural Si contains about 5% 29 Si or less through isotopic 
purification. Isotopic purification or nuclear polarization will 
suppress spectral diffusion in Si. 



quantum computation, where the localized electron spins 
are spatially well separated, a direct magnetic dipolar in- 
teraction between the electrons themselves is not an im- 
portant dephasing mechanism^ Therefore, the interac- 
tion between the electron spin and the nuclear spin bath 
is the important decoherence mechanism at low temper- 
atures and for localized electron spins. Now we restrict 
ourselves to a situation in the presence of an external 
magnetic field (which is the situation of interest to us) 
and consider the spin decoherence channels for the local- 
ized electron spin interacting with the lattice nuclear spin 
bath. Since the gyromagnetic ratios (and, hence, the Zee- 
man energies) for the electron spin and the nuclear spins 
are typically a factor of 2000 different (the electron Zee- 
man energy being larger), hyperfine-induced direct spin- 
flip transitions between electron and nuclear spins would 
be impossible (except as virtual transitions as will be 
discussed in Sec. Ill Ap at low temperature since phonons 
would be required for energy conservation. This leaves 
the indirect SD mechanism as the most effective electron 
spin decoherence mechanism at low temperatures and fi- 
nite magnetic fields. The SD process is associated with 
the dephasing of the electron spin resonance due to the 
temporally fluctuating nuclear magnetic field at the lo- 
calized electron site. These temporal fluctuations cause 
the electron spin resonance frequency to diffuse in the fre- 
quency space, hence the name spectral diffusion. These 
fluctuations result from the dynamics of the nuclear spin 
bath due to dipolar interactions with each other along 
with their hyperfine interactions with the qubit. This 



scenario is illustrated in Fig. [TJ Nuclear spins in the spin 
bath flip-flop due to their mutual dipolar coupling (since 
the typical experimental temperature scale ^100 mK is 
essentially an infinite temperature scale for the nuclear 
spins with nano-Kelvin scale coupling), and this leads to 
a temporally random magnetic-field fluctuation on the 
central spin, i.e., the electron. 

A. Interactions 

In the most general form, the SD Hamiltonian for the 
central spin decoherence problem may be written (in h — 
1 unit) 

n = n e s z + £ /„ ■ A n • I (i) 

n 

"F ^ ^ LUnlnz "F ^ ^ In ' B nm ' Ijni 
n n<m 

where S and /„ are vectors of spin operators for the elec- 
tron (central spin) and nucleus n respectively, VL e and 
uo n are their respective Zeeman energies (with an ex- 
ternal magnetic field applied in the z direction), and 
A n and B nm are tensors describing electronuclear and 
internuclear interactions respectively. Isotropic Fermi- 
contact hyperfine (HF) interactions typically dominate 
A n (i.e., A n — A n t) although anisotropic FIF inter- 
actions, due to dipolar contributions, may also be im- 
portant. Internuclear dipolar interactions often domi- 
nate B nm , though other local interactions between nuclei 
such as indirect exchange interaction o 26 i 27 ' 28 may also be 
significant £ Typical energy scales are shown in Table U 
for convenience. 



TABLE I: Interactions and estimated energy scales. (A simi- 
lar table appears in Ref. Q.) 



Interaction 


Symbol 


Scale (h = 1) 


Scale (fee = 1) 


Zeeman (electron) 




10 11 s" 1 a 


1 K a 


Zeeman (nucleus) 




10 8 s" 1 a 


1 mK a 


Contact HF 


A n 


10 6 s" 1 6 


10 fiK b 


Cumulative HF 




< 10 11 s" 1 c 


<1K C 


Dipolar 




10 2 s" 1 


1 nK 


Indirect exchange 


R Ex 


10 2 s" 1 


1 nK 


HF mediated 


A nm 


10 s" 1 


1Q- 1 nK 



a With a magnetic field of about 1 T. 

b In III-V compound quantum dots with N ~ 10 5 nuclei and 
also in Si:P donors. 

c In III-V compound quantum dots, A ~ 10 1 s . In Si:P 
donors, A ~ / x 10 10 s _1 , where / is the fraction of 29 Si. 



The HF energies are typically many orders of magni- 
tude larger than internuclear dipolar energies: ||A„|| 3> 
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II B nm || . By ignoring the B nm term for a moment, de- 
coherence may occur as a result of real or virtual elec- 
tronuclear flip-flops via the HF interaction; 7 -^ Such a 
process may be suppressed, however, by increasing the 
applied magnetic field due to a conservation of Zccman 
energies. These Zeeman energies are given by 51 e = 7s-B 
and w„ = — 7 n -B, where B is the applied magnetic-field 
strength, and 75 and "f n are gyromagnetic ratios of the 
electron and nucleus n, respectively (with 7„ defined in 
an opposite sense as 75). Typically, 75 ~ 10 3 7„ so that 
the nuclear Zeeman energy is negligible relative to the 
electron's Zeeman energy and the electron must over- 
come its Zeeman energy barrier in order to flip. In the 
limit of O e 3> A, with A = J2 n A n , HF-induced elec- 
tronuclear flip-flops are effectively suppressed. When we 
reintroduce B nm , however, decoherence still occurs and 
is well defined in the f2 e — > 00 limit because the electron 
spin will dephase as a result of nuclear field fluctuations 
induced by inter-nuclear interactions. This decoherence 
is spectral diffusion. 

In the limit of a large applied field (formally we will 
say tt e — > 00), electron flips are completely suppressed. 
In this limit, the effective Hamiltonian becomes TL ~ 



J2n W n^n2 + En 



n<m I n - B nrn ■ I m , where 

A n is the z column vector of A n . We are free to drop 
Vt e S z for any dynamical considerations now because it is a 
conserved energy in this limit. Effects due to anisotropic 
HF interactions may be treated independently of SD with 
a trivial disregard to internuclear interactions^ For our 
purposes, therefore, we will treat only the z component of 
A n from the Fermi-contact HF interaction, which leaves 
us with 

T~L ~ ^ A n I nz S z + ^ ' uj n I nz + ^ I n ■ B nm ■ 7 m ,(2) 

n n n < rn 

A n = y 7S 7nft|*(fin)| 2 . (3) 

Equation ([3]) gives Fermi-contact HF coupling constants 
that are proportional to the probability of the electron 
being at the nuclear site (the square of its wavefunction) . 

How large of a magnetic field must we apply for Eq. © 
to be a valid effective Hamiltonian? To be explicit, 
one may use quasidegenerate perturbation theory 3 "^ to 
systematically transform Eq. ([1]) into the block form 
TL s» E ± |±)7i±(±|, where {|+),|— )} are the electron 
spin states. This transformation will be convergent if 
il e A, where A = J2 n is the maximum cumulative 
HF field. In III-V compound quantum dots, fi e ~ A with 
an applied field of 1 T, calling this approach into ques- 
tion; however, the Q e 3> A condition is overly strict be- 
cause this maximum HF field is never reached in a bath 
that is not fully polarized. It is more relevant to con- 
sider how higher orders of this transformation [Eq. @ 
represents the zeroth order] compare with the internu- 
clear Bnm coupling term to generate decoherence. In the 
next order of the transformation, HF-mediated interac- 
tions between nonlocal nuclei emerge. This interaction, 



wellknow n 31 ! 32 as the Ruderman-Kittel-Kasenya-Yosida 
interaction, results from virtual electron spin flips and is 
suppressed with a large applied magnetic field. This next 
order contribution to the Hamiltonian is given byS. 



TL 



HFM 



E 

n^rn 



■A-rim In _|_ I-m — $z 



(4) 



where A nm — A n A m /VL e . In applying this transforma- 
tion, we must slightly rotate the basis states; this re- 
sults in a "visibility" loss 33 of coherence estimated as^ 
E n (A n /Qe) 2 and is certainly small for a large bath in 
which £! e > A. In this paper, we will restrict ourselves to 
the £l e — > oo limit and merely mention how HF-mediated 
or other higher order interactions may play a role. 

Direct interactions between nuclei are often dominated 
by their magnetic dipoles with the form 



H 



"r 3 



R 5 



. (5) 

The dipolar interaction between nuclear spins in semicon- 
ductors has a typical strength of 7i^ TO ~ 10 2 s _1 , which 
is much smaller than typical nuclear Zeeman energies of 
about 10 8 s _1 in an applied field of 1 T. Therefore, 
energy conservation arguments allow us to neglect any 
term that changes the total Zeeman energy of the nu- 
clei. This leaves the following Zeeman energy-conserving 
secula r 31 ' 34 part of the interaction: 



i 



E B " 



n < rn n^m 



11-3 cos 2 1 



(7) 



where R n m is the length of the vector joining nucleus n 
and nucleus m, and 9 nm is the angle of this vector rela- 
tive to the z magnetic-field direction. The 8 



7n,7r. 



denotes 

that the flip-flop interaction between nuclei with differ- 
ent gyromagnetic ratios should be excluded because of 
Zeeman energy conservation in the same way that the 
nonsecular part of the dipolar interaction is suppressed. 
This occurs, for example, in GaAs because the two iso- 
topes of Ga and the isotope of As that are present have 
significantly different gyromagnetic ratios. This secular 
interaction corresponds to a B n m matrix for Eq. ©, with 
B nm , Bnm, and —2B nm along the diagonal, respectively, 
in the x-y-z spin basis for a pair of nuclei having the same 
gyromagnetic ratios. (The I nz I m z term plays an insignif- 
icant role, which is why we use the same B nm symbol, in 
different fonts, for the scalar and tensor). 

In addition to the dipolar interactions between nuclei, 
an indirect nuclear-spin exchange ; 26 ' 27 ' 28 ' 32 ' 35 which is 
mediated by virtual electron-hole pairs, may also have a 
significant quantitative impact on SD in III-V materials. 9 
This interaction takes the form 



a/ax R 1 ^ J 

^nm nm n 



(8) 
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The corresponding B nm is B„* = The leading 

contribution to this pseudoexchange interaction for near- 
est neighbors may be expressed a o 26 ' 35 



dEx __ 



Ex„,Ex 



MO In % 



ao 



i? 3 E 



(9) 



where 7^ x is the effective gyromagnetic ratio determined 
by a renormalization of the electron charge density^ This 
interaction has been experimentally studied many years 
ago . 26 ' 27 i 28 In GaAs quantum dots, these interactions can 
be comparable to the direct dipolar interactions. There 
may be other local interactions between nuclei in the 
bath, such as the indirect pseudodipolar interaction^ 2 , or 
internuclear quadrapole interaction, but the dipolar and 
indirect exchange interactions alone account for nuclear 
magnetic resonance line shapes^ i 26 ' 27 i 28 In any case, all 
such local interactions may easily be included in our for- 
malism. 

To summarize and put everything in a convenient gen- 
eral form, we approximate our Hamiltonian as Tt « 
£±|±>W±<±| where 



(10) 



TL e b is the electron dependent part that plays the role 
of coupling the electron spin to the bath, and Tibb in- 
cludes secular (preserving nuclear Zeeman energy) bath- 
bath, i.e., internuclear, interactions such as the secular 
dipolar interaction [Eq. ([6])] and the exchange interac- 
tion [Eq. ([5])]. Zeeman energies are omitted from this 
Hamiltonian because they are preserved by all included 
interactions and, thus, not relevant to the dynamics. 
In the fL — > oo limit, 



H 



eb ~ TT 



1 



2 ^ 



A T 



(11) 



with a factor of 1/2 from the magnitude of the electron 
spin, and A n given by Eq. ([3]) with a proportionality to 
the square of the electron wavefunction at site n. To vali- 
date this approximate form, one may consider the effects 
of higher order interactions of the canonical transforma- 
tion, such as the HF-mediated interactions [Eq. Q] that 
contribute to H e b (due to its S z dependence). These 
higher order interactions introduce additional decoher- 
ence mechanisms that will shorten coherence times at 
lower magnetic fields. The main consideration of this 
paper is the decoherence that may not be removed by 
simply increasing the magnetic-field strength and is ap- 
plicable in the limit of a large applied field at the point 
where decoherence is insensitive to the strength of the 
applied field. 



Since we consider a dephasing-only Hamiltonian of the 
form H ~ J2± l = ' = )^± we need only deal with dephas- 
ing decoherence. Dephasing decoherence involves only 
the transverse component of the electron spin. For a 
given experiment, we define the up and down evolution 
operators, U ± , as evolution operators for the bath given 
an initially up or down electron spin. If we just have free 



evolution for a time t, then U± = Uq = exp 



C-iH±t 



In 



general, we can consider an experiment with a sequence 
of 7r pulses that flip the electron spin between periods of 
free evolution time t,-, so that 



U± = ... exp (— iH T T2^ exp (— iH±Ti\ . 



(12) 



The transverse component of the expectation value of the 
electron spin will then decay, in magnitude, by a factor 



of 



(W) 



where W 



an appropriately weighted average over the bath states, 
and ||...|| takes the magnitude of the resulting complex 
number. The coherence decay is thus characterized by 

(W) . In an echo experiment, one initializes an ensem- 
ble of pure electron spin state in some transverse direc- 
tion, applies a sequence of it pulses designed to refocus 

the spins, and observes an echo signal, ve = (W) , at 

the end of the experiment. 

Given an arbitrary initial bath density matrix written 
in the form pb = Ylj Pj\Bj)(&j\> we m ay average W over 
bath states with 



denotes 



(W) = J2Pj(Bj\W\B j ) 



(13) 



where each \Bj) is a different bath state. By referring 
to the energy scale estimates of Table HI temperatures 
on the mK — K scale are justifiably treated as infinite 
with respect to the nK scale intrabath interactions. The 
remaining Zeeman and Fermi-contact HF interactions in 
our approximate Hamiltonian quantize the nuclear spins 
in the z direction (involving only I nz nuclear spin op- 
erators) and, thus, we can approximate the initial bath 
density matrix as a mixed state composed of uncorrelated 
pure nuclear spin states in this basis so that 



\B, 



n 

n 



Prim | In j ^) 



(14) 



where I n is the magnitude of the nth nuclear spin and 
\I n ,m) n represents the state of the nth nuclear spin 
with a z projection of m. The cluster approximation in 
Sec. IIIIAl will make use of the assumption that the bath is 
initially uncorrelated (at least, approximately). Section 
IIVI will use the z quantization as a further convenience. 



B. Decoherence 



C. Dynamical decoupling pulse sequences 



We characterize decoherence as the expectation value 
of the electron spin over the evolution of the experiment. 



If we simply let the system freely evolve, the decay of 
the electron spin expectation value strongly depends on 
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the type of averaging we perform over bath states. If 
we consider an ensemble of electron spins, each with its 
own bath, then we will see rapid electron spin dephas- 
ing simply due to the distribution of HF nuclear fields, 
'^2 n A n I nz /2. This is known as inhomogeneous broad- 
ening because it broadens the electron spin precessional 
frequency due to inhomogeneity of the effective magnetic 
field. This, however, is an artifact of the ensemble or our 
lack of knowledge of the effective nuclear field of a static 
bath and not true decoherence. If we consider a single 
electron spin with a known nuclear field, or a homoge- 
neous ensemble which may be obtained through mode 
locking^ for example, we arrive at a dynamical decoher- 
ence known as FID. In our dephasing model, FID can 
only arise from interactions among bath elements, such 
as local dipolar or nonlocal HF-mediated interactions. 

Traditionally— nonhomogeneously broadened coher- 
ence is measured from Hahn spin echoes. The Hahn 
echo sequence simply involves a single n rotation mid- 
way through the evolution such that U±(t) — U^{t) = 
U^(t)Uq(t), with r = t/2. We denote this sequence 
with t — > 7r — > t\ free evolution for a time r, then a 7r 
pulse, then free evolution for a time r again, with the 
arrows indicating the sequence ordered in time. This 
sequence will reverse the effect of any inhomogeneous 
static field. What remains is SD induced by a dynami- 
cal nuclear bath. It is important to note, however, that 
the effects of the Hahn echo go beyond the elimination 
of inhomogeneous broadening. The Hahn echo is also 
a DD sequence,— in which the first order of the Mag- 
nus expansion^ is removed by the fact that the time- 
averaged Hamiltonian, proportional to TL + + 7i_ oc 7if,, 
decouples the qubit from the bath. For this reason, the 
Hahn echo does not have the same effect as homoge- 
neous (or single qubit) free induction decay— In partic- 
ular, the Hahn echo removes the lowest-order effects of 
HF-mediated interactions. 9 This is because HF-mediated 
interactions, having an S z factor [Eq. Q], belong to H. e b, 
and if we consider only HF-mediated intrabath interac- 
tions, then fj^ = cxp (^kiH^T^j exp (jfiH e bT^ = 1- 

Given a DD sequence, such as the Hahn echo for a 
dephasing system, coherence over a given net amount of 
time may be increased through a rapid repetition of the 
basic sequence. This strategy, known as bang-bang in the 
quantum information community," gives coherence en- 
hancement at the cost of more frequent applications of 7r 
pulses. Pulses must be applied more frequently because 
errors due to higher order terms of the Magnus expan- 
sion pile up over the course of the sequence. A better 
strategy is to use recursion, rather than repetition, to 
generate concatenated sequences^ Such CDD, with the 
Hahn echo as the base sequence, has been showniS to be 
effective for the SD problem. In fact, with each concate- 
nation, we demonstrated, in Ref. Il5l . coherence enhance- 
ment with an increase in the time between pulses (let 
alone, the net sequence time). 

With I > levels of concatenation, our CDD pulse 



sequence is recursively defined by-2 9 - p ; := — > n — > 
Pi_i — > 7T, with p := r. CDD with I = is simply free 
evolution. At I = 1, we have p x := r — > 7r — > r — > tt, 
which is simply the Hahn echo (with an extra ir pulse 
at the end to bring the electron spin back to its original 
phase apart from the decoherence). With each concate- 
nation, we do to the previous sequence what the Hahn 
echo does to free evolution and, in this way, we obtain 
improved DD. This sequence may be simplified by noting 
that two 7r pulses in sequence do nothing. Thus, 



Pi 




, odd I 
x , even I 



(15) 



(again, arrows indicate sequences ordered in time) and 
the up and down evolution operators at level I have the 
recursive form ofi^ 



(16) 



Recently, a series of DD sequences was discovered by 
Uhrig 1 - to be optimal in the number of pulses for the 
spin-boson model. The n-pulse sequence in this series 
may be defined by 



7r(j - 1) 

cos I I — COS 

n + 1 



(17) 



for 1 < j < n + 1. Unlike CDD, Uhrig DD (UDD) re- 
quires only a linear (rather than exponential) overhead 
in the number of pulses for each order of coherence en- 
hancement. Furthermore, UDD was showni& to kill off 
successive orders in a time expansion in a completely 
model-independent manner. The UDD sequence has a 
strong advantage over CDD in its linear versus exponen- 
tial scaling of the number of pulses; however, it is effective 
only when a time expansion is convergent, while CDD is 
also effective in the intrabath perturbation for SD (this 
is important since B nm -C A n ). 

In Sec. HVi we will discuss the wavefunction depen- 
dence in the short-time approximation. It would be nat- 
ural to discuss this in the context of the UDD series since 
UDD is effective in this short-time limit. However, the 
methods of Sec. II V CI work for the CDD series but, un- 
fortunately, do not work for the UDD series. In this 
work, we therefore focus attention on FID and the CDD 
series. Note, however, that the UDD and CDD series 
are the same for levels zero (FID), 1 (Hahn 1 ^), and 2 



III. CLUSTER METHOD 

In this section, we review our cluster expansion 
method&ii for solving the SD problem in a particularly 
simple and illuminating form. Section [ill Al gives the ba- 
sic cluster expansion result, the cluster approximation, in 
which we equate (W) [see Sec. Ill B| to the exponentiation 
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of single-cluster contributions. This is useful because per- 
turbation theory cannot be directly applied to (W), but it 
can be applied to its single-cluster contributions (it is the 
multicluster contributions that are particularly problem- 
atic for perturbation theories due to the large number of 
possibilities in combing different clusters). We consider 
aspects of such perturbation expansions in Sec. lIII Bl that 
will be relevant to Sec. ITVl 



A. Cluster approximation 

The cluster expansion is based on the fact that our H± 
Hamiltonians couple nuclei via relatively weak pairwise 
interactions in a large bath with N nuclei (any n-way in- 
teractions could justify such an approximation as long as 
n <C N). The W operator, which is the product of evo- 
lution operators arising from Tl±, can, in principle, be 
expanded into a sum of products of the Hamiltonian in- 
teraction elements. The interaction elements of each such 
term will uniquely determine a set of clusters of nuclei; 
each cluster, with respect to this term, may involve inter- 
actions among itself but not among any other cluster, and 
no cluster may be divided into further subclusters. For 
example, a term with BipIi+h—Ba^Ia+Ii— forms two 
clusters, {1,2} and {3,4}, while Bi a h+h-B 2 .s,h+h- 
forms a single cluster of {1, 2, 3}. The term "cluster" im- 
plies proximity among the member nuclei as applicable 
to local dipolar interactions; however, we may also treat 
nonlocal HF-mediated interactions using the term clus- 
ter in a more general sense as a set of nuclei that are 
interconnected by interactions under consideration. 

If one considers a perturbative expansion of (W) with 
respect to the pairwise interactions, one immediately 
faces the problem that the number of terms of W scales 
in powers of ./V with successive inclusion of the pairwise 
interactions, destroying any hope of convergence when N 
is large. To resolve this problem, let us first partition W 
according to the number of clusters involved in each term 
such that 



JV 



w 



p=0 



(18) 



where we define as the sum of the terms from W that 
involve p independent clusters. Note that = 1. By 
considering only local interactions (e.g., dipolar), then it 
is apparent that a perturbative expansion of (pyM) with 
respect to the pairwise interactions does not suffer from 
the adverse N scaling suffered by (W) because, when in- 
teractions are local, there are O (N) clusters of any size. 
Even with nonlocal interactions, the perturbative expan- 
sion of (W - ! 1 !) is a generally significantly better controlled 
expansion for large N than that of (W). 

Assume that, arising from a perturbative expansion 
with respect to the pairwise interactions, is well 

approximated when only including contributions due to 



clusters of some maximum size that is much less than 
N. Along with our assumption that the bath nuclei are 
initially uncorrelated, Appendix lAl shows that (W^) w 
(W^) k /k\ [Eq. (|A5[) ]. This approximation breaks down 
as k becomes significant relative to N; however, for a 
large enough bath where the previous assumptions are 
met, 



(W) « exp 



(<^W>) 



and 



v E 



(W) 



exp { Re 



.{<#M>}). 



(19) 



(20) 



A formalized cluster expansion, with a discussion on con- 
vergence tests, is presented in Ref. [ll|. The cluster ap- 
proximation presented in this subsection gives an equiv- 
alent result with a simpler derivation. 



B. Single-cluster perturbation 



Given the result of Eq. pOjl , we have reduced the prob- 
lem of SD decoherence to that of perturbatively treating 

Re |. That is, we wish to consider a perturba- 

tive expansion of W where we neglect terms involving 
multiple clusters. To be consistent with the cluster ap- 
proximation, such a perturbation should be directly or 
indirectly tied to cluster size so that large clusters may be 
neglected. To this effect, wc may perturbatively treat the 
pairwise interactions (the intrabath perturbation) or we 
may consider a time expansion which, in a sense, pertur- 
batively treats all of the interactions (e.g., HF, dipolar, 
and HF mediated); in either case, larger clusters require 
more interaction factors and thereby increase the order 
of the perturbation. 

In this section, we will consider general perturbative 
properties that apply to both the intrabath and time 
perturbation expansions. In Sec. IIV1 we will specifically 
consider the time perturbation and see how it may be 
used in the formulation of a convenient continuum ap- 
proximation. In a perturbation whose order increases 
with increasing cluster size, the lowest order of (W) — 1 
is equivalent to the same order of (W^ 1 ') because terms 
of W with multiple clusters are automatically higher or- 
der terms (products of lower order terms). To the lowest 
order, then, and in the context of CDD with I levels of 
concatenation, we will consider (W). 

By noting that C/j are unitary operators such that 

t___ 

we may write 



U, 



Re 



{<*>} 



u; 



u. 



2X AtA, 



U- ) (21a) 
(21b) 
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where we define A ; = U+ - Uf. Thus, ( AjA^ gives a 

measure of the decoherence. By applying the recursive 
definitions for the fjf^ evolution operators [Eq. jT6])] , 



A 



Uf_ v U+ 1 = Uf-v&l-i , (22) 



since U,^ commutes with itself. 



Let us consider a perturbation with a smallness param- 



eter A in which U, 



O (A) for all I > Iq for some 



Iq. Two such perturbations are the time expansion with 
A = r, Iq = (since no evolution occurs in the r — > 
limit) and intrabath perturbation with X — e, l = 1 
(since there is a perfect spin echo refocusing in the e — > 
limit). Because the identity commutes with anything, it 
is easy to see from Eq. (l22l that A; = (A) x A;_i for 
all I > l ; this proves that we get successive cancellations 
of the low-order perturbation (r or e) with each concate- 
nation of the sequenced The lowest order result is given 
by 



d - 

dX U ^ 



A=0 



, V I > l . (23) 



Conveniently, for all I > Iq, 



dA 



so that Eq. 



A=0 



~ d\ Ul - 1 



(24) 



= 2 l -> 
becomes 



4r (uf + ur ) 12 

dX\ l ° l ° 1 ' 



A=0 



A=0 



l-lo-l 



2«-j 



. A;_i 



A=0 



, l>lo 



A=0 



, I = Iq 

(25) 



Note that in the Z = Zq + 1 case, Eq. (|22|) yields 



A Jo+1 « A^ 



d - _ d » i 
dX l °'dX l ° 



(26) 



A=0 



which is equivalent to the corresponding case in Eq. (I25p 
recalling that any operator commutes with itself. 



IV. WAVEFUNCTION DEPENDENCE IN THE 
SHORT TIME LIMIT 

In this section, we use the formalism developed in 
Sec. IIII1 with the cluster approximation and general per- 
turbation formulation of Sec. IIII B| to derive results ap- 
plicable in a short-time limit and use these results to 
formulate a continuum approximation useful for under- 
standing the dependence of spectral diffusion on the 
shape of the electron wavefunction. In Sec. IIV A[ we ap- 
ply the general results of Sec. IIII Bl to obtain the lowest 
time perturbation results of Re{(M / [ 1 ')}. This is done 



for the cases of free induction decay and concatenated 
echoes. Section IIV Bl shows how the nuclear dependent 
and electron dependent parts of this lowest-order time 
perturbation solution may be separated in a way that 
allows us to generically treat the bath for any electron 
wavefunction. Scction llV Cl takes this one step further by 
treating the bath as a continuum so that we may obtain 
results via integration for any given electron wavefunc- 
tion. In Sec. IIV Dl we discuss the circumstances in which 
the time expansion may or may not be applicable. 



A. Time perturbation 



Section IIII Bl considered the perturbations of concate- 
nated Hahn sequences in a general sense in the context 
of the cluster approximation of Sec. IIII Al Now we are 
specifically interested in the time expansion. We first 
address this in the case of FID and then treat the Hahn 
sequence and its concatenations. 



1. Free induction decay 

The I = result of A = U + — [/_, without any pulses, 
gives A = -2iH eb T. Where H eb = J2n A ^^ lz /"2, this 
result for Ao is simply due to inhomogeneous broaden- 
ing. In the case of free induction decay, we are not con- 
cerned with inhomogeneous broadening and would like to 
obtain the SD decoherence of a single electron spin (or a 
mode- locked ensemble^). We may do this by entering the 
rotating frame of reference for the electron precessing in 
the nuclear field. This may be done, effectively, by mak- 
ing the transformations, U^t) — > U^^exp (^tiii^t 



for free induction decay, we then obtain 



An = 



7~Leb,7~Lb t . 



(27) 



The free induction decay is then given by Eq. (12T))) via 
Eq. (|21b|) with this form of A . 

2. Concatenated dynamical decoupling 

We now consider the Hahn echo sequence and its con- 
catenations: I > 1. There is no need to go into a rotating 
frame as we did for FID because these sequences auto- 
matically reverse the effects of inhomogeneous broaden- 
ing. We will apply Eq. P5|) by using A = r. In the limit 
of r — > 0, no evolution can occur, so it is apparent that 



Uq = 1 and, thus 

T = 

the Hahn echo [Eq. 



i = in this context. For I 
becomes 

HebiHb 



(28) 



Note that there is a simple relationship between Eqs. 
and (|27p for the Hahn echo and FID, respectively. The 
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relationship is not so simple when we move away from the 
Q e — ► oo limit and consider HF-mcdiated interactions. 
In that situation, discussed in Ref. 0, the lowest order 
contribution to FID in the time expansion will come from 
HF-mediated interactions but this lowest order effect will 
be cancelled in the Hahn echo. 

To consider concatenations of the echo, we simply ap- 
ply the recursion of Eq. (|25p to obtain 



_ 2 (; 2 -;+2)/2 



T~(-eb, T~(-b 



(ir) 



1+1 



0(t 1+2 ), VZ>0, 



(29) 



with / nested commutations abbreviated by ...'s. As a 
result of these nested commutations and as observed in 
Ref. [HI, each concatenation introduces larger cluster sizes 
to the lowest-order expression (i.e., each time we com- 
mute with H.b, we may introduce an additional nuclear 
site to any term of this operator). 



B. Pair "correlations" 

A reasonable assumption for many solid-state spin 
baths is that the bath Hamiltonian Tib, which excludes 
qubit-bath interactions, is homogeneous. That is, sites 
that are equivalent in terms of the Bravais lattice are 
equivalent with regard to bath interactions. A notable 
exception to this is where isotopes in the lattice are in- 
terchangeable; for example, three different isotopes of Si 
may occupy any lattice site in Si, and two different iso- 
topes of Ga may occupy the Ga sublattice in GaAs. How- 
ever, if we simply want to know the decoherence that re- 
sults from averaging different types of isotopic configura- 
tions, then we may regard the bath (apart from the qubit 
interactions) as homogeneous and use isotopic probabil- 
ities in expressions for Tib- Then the only inhomogene- 
ity is in the interactions with the qubit, 7i e f>- We can 
then factor out this inhomogeneous part and compute 
the rest in a way that is independent of the qubit in- 
teractions. This will be convenient, for example, when 
analyzing a quantum dot in which the wavefunction of 
the electron (whose spin represents the qubit) can take 
on many shapes and sizes. 

By referring to Eqs. (j2"Tj) and (j2"9"|) , we can make the fol- 
lowing factorization of the homogeneous and nonhomo- 

geneous parts of ( A J A; \ [determining SD via Eqs. (|21bp 



and PU| ] 



-(AfA* 



EA* A f(0 t p 
+ O (rf+ 2 ) , 



(30) 



For I = (FID), 



f(0) 



F(o) = 



F. 



(<>) 



rnz I ? 



Inz , Hb 



and for I > 0, 

/W = (_)('+i)2( /2 -' +1 ) 



pin 



Inz, Hb 



(32) 
(33) 

(34) 
(35) 



where the ...'s again denote I nested commutations. 

Since we assume the high field limit where secular cou- 
pling Hb preserves nuclear polarization, Hb,J2 n ^ nz = 
so that J2 m fn,m = J2 m fm,n = for any n. Then, 

/n,n = - J2mjtn (/n,m + fm,n^ /2. By using this fact, 

we may rewrite Eq. (|30|) in terms of the differences of the 
HF constants; after all, if A n is the same for all nuclei, 
there is no nuclear induced spectral diffusion in the high 
field limit (nuclear flip-flops would have no effect on the 
electron) . We will assume that the HF constants are real, 
A n = A* n , as is the case for the Fermi-contact interaction 
[Eq. ©]. Then, 



^ ] A n A m f, 



(0 

n.m 



l^2A n A m ^l + fg] n ) (36) 
+ Vi 2 f (,) 

n 

-- V (A - A ) 2 ( f (l) + f (i) ^ 

i / j X^-n s±m) \Jn,m • J m,n I • 



Thus, in the short-time limit of Eq. (f30|) . 



2(a;a ; 



-i^(A„-A m ) 2 Re{/« m p. (37) 



The homogeneous part is represented by /n,L, and by 
exploiting this homogeneity, we note that this function is 
equivalent when we shift by any Bravais lattice vector, R: 
/« (r n , r m ) ee fi% = fW (f n - R } f m -R). We may then 

relate any r n — R to one of the basis sites of the Bravais 
lattice, so then, with b representing the corresponding 
basis site of r„, we may write 



where 



/ra,m — fb( r n r m)- 



C. Continuum approximation 



(38) 



p(0 



4 ,1 = 
21 + 2 , I > 



^32) The pair correlation formulation above is particularly 

convenient in the context of a continuum approximation 
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for HF coupling constants. From the Fermi-contact HF 
interaction [Eq. (J3J], 



A =^ 77 

3 \ a 6 



P(f n ), 



(39) 



where d n is the charge density for the isotope at site n, 
V u is the volume of the unit cell of the Bravais lattice, a is 
the lattice constant, and P(r) oc ||vl/(r)|| 2 is the electron's 
probability density normalized such that J d f f P{r)/a — 
1. 

Let £ characterize the correlation length scale from 
Eq. (O. Then, if \diP\ / |<9 2 P| « i, where d t = d/dxi 
and (m — (d/dxi) (d/dxj), for a given n and m pair in 
Eq. I30l we may use 

P(f n ) - P(r m ) » (r„ - r m ) • VP (f m ) . (40) 

Furthermore, by using a continuum approximation where 
we replace one of the summations with an integral, 
Eq. I|37p then becomes 



1 
2 

where 



Y 



£3 



Cf=n c 



'-j el i fid I 



o 5 " 



r 3 -VP(f) 



(41) 



1 2 



(42) 



and n c is the number of sites in the conventional cell 
(of volume a 3 ) and (...)/ averages over different isotopes. 
Note that n c is equal to the number of basis sites multi- 
plied by a 3 /V u . 

By rearranging Eq. (|41|). 



AjAi) « r" / ^ [vP(r)] T M^ [vP(r) 



(43) 



where x 1 denotes the transpose of any vector x, [fj] [rj] T 
is an outer product, and M"' is a matrix. It is im- 
portant to note that M^- 1 is independent of the electron 
wavefunction (or its probability density); these are con- 
stants that are predetermined for a particular lattice and 
applied magnetic-field direction. The wavefunction de- 
pendence is entirely of the form j d 3 r[diP(r)] [djP(r]]. 
Being symmetric, we may diagonalize to the form 
M« = u t ufm ( p so that 



-(A| A, 



rf 3 



/ — u t -VP(r) 



(44) 



Details of how we actually computed M^> for various 
systems are given in Appendix [Bj 



Putting this in yet another form, 

p(0 



i(At A; ) = [r/4 l f\o(r^), (45) 

r -> i 2 
1 / d 3 r itj • VP(r) /a 



.(0 



.(0 



where p(l) is defined by Eq. (|31|1 and the = 

-l/(2/+2) 

have units of time. For a given con- 
catenation level I, the echo signal [Eq. (|20|) ] is approxi- 
mately 



(0 / 2 



ve ~ exp 



p(0 



(47) 



in the limits of a strong applied magnetic field and in 
the short-time approximation (i.e., extrapolated from the 
short-time behavior which may or may not be valid at 

T = T ). 

We will now treat, specifically, the case of a quantum 
well with thickness zq and Fock-Darwin radius ro (re- 
sulting from a combination of parabolic confinement and 
confinement due to the magnetic field) £ In this case, the 
wavefunction is sinusoidal in the z direction and has a 
Gaussian form in the lateral direction. The probability 
density is of the form 

x 2 + y 2 " 

P{x,y,z) cx exp ( * — I \ 



cos 



(ji^e(V2-M)- (48) 



Let us consider the case where the z vector is an eigen- 
vector of (e.g., when the problem, with the applied 
magnetic field direction, is symmetric about the z axis). 
By recalling that P{r) should be normalized such that 
fd 3 rP(f)/a 3 = 1, 



/ 



d 3 f 



~dP(r)' 


2 = S' ef 


~dP(r)~ 


dx 




dy 




3a 6 






Anzor^ 




~dP{r)~ 


2 ira 6 




dz 


4 r o' 





d 6 r 

We may simply plug these into Eq 

-p(0 3(a/r ) 



(0 _ j' ttQ/zq) 
_i N n 



(0 



2nN n , 



(49) 



(50) 
as follows: 



(51) 

defining N a = V/a 3 as the number of a x a x a lattice 
cubes in a quantum dot of volume V — z^r 2 . When 
t"o 3> Zo> for example, we have 



.(') 



,(0 



A^o/a) 2 



(52) 
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Table |TT] shows computed values of fi • ' for the GaAs 
or InAs lattice with < I < 3 (free induction decay, 
Hahn echo, and up to three levels of concatenation). Sub- 
scripts of n in this table indicate corresponding eigenvec- 
tor, Hi, lattice directions. GaAs and InAs both have 
a zinc-blende structure with the 75 As atoms on one of 
the two fee lattices and respective lattice constants of 
5.65 and 6.06 A. — The natural abundance of Ga isotopes 
60.4% 69 Ga and 30.2% 71 Ga, and the natural abundance 
of In isotopes are 4.3% 113 In and 95.7% 115 In42 The gy- 
romagnetic ratios are 77 = (4.60,8.18,6.44,5.90,5.88) x 
10 3 (s G)" 1 for 75 As, 71 Ga, and 69 Ga, 113 In, and 115 In 
respectively.^ 3 * We have also used the following respective 
charge densities d(I) = (9.8, 5.8, 5.8, 2.3, 2.3) x 10 25 cm" 3 ; 
these charge densities were estimated in Ref. for GaAs 
and estimated in Ref. using the technique in Ref. 
for InAs. The Ga and As nuclei have spin magnitudes 
of / = 3/2 and the In nuclei have spin magnitudes of 
I = 9/2, which we account for properly. The table shows 
the results for an applied magnetic field in the [001] or 
[110] lattice directions. For the bath-only Hamiltonian 
TLb 1 we have included the secular dipolar interaction 
[Eq. ©] and, for GaAs, indirect exchange interaction 
[Eq. ®] witb22 B% = ~ lnlm h(V^6 A/2iZ* m ). The 
table shows GaAs results when we include or exclude 
indirect exchange interactions; the remainder of the ta- 
ble only considers spectral diffusion induced by dipolar 
interactions among bath nuclei. Because of the near de- 
generacy of the In gyromagnetic rations and the natural 
predominant abundance of 115 In, we simply show results 
for 100% 115 In; this should also yield the lower bound of 
dipolar-induced spectral diffusion decoherence times for 
any InAs/GaAs mixture (since In induces the strongest 
decoherence due to its large 7 = 9/2 spin). Mixing a 
little Ga into InAs or even mixing a little In into GaAs 
will increase decoherence times as a result of the reduced 
probability for any given cluster to be of the same nuclear 
isotope (only these clusters can contribute in the high 
magnetic field limit). Contributions from the As nuclei 
to decoherence is unavoidable in any such mixture; we 
thus give As-only results in the table as an upper bound 
for decoherence times in InAs/GaAs mixtures. This As 
contribution, however, will vary depending upon lattice 
constant; we give the full range in which slightly longer 
decoherence times result from using the slightly larger 
InAs lattice constant and vice- versa for the GaAs lattice 
constant. 

It is important to note that the decay time for the 
overall CDD pulse sequence at level I is to = 2 1 tq. Values 
of r represent the time between pulses rather than the 
overall time. Thus, the fact that the values of /Xj [related 
to To via Eq. (|4"6"|) ] increase in Table |TT] as / increases 
beyond I = 1 yields extra coherence time enhancement 
beyond the 2 l extension of to- That is, as noted in Ref. [HI, 
not only does concatenation increase the net coherence 
time, but it also decreases the frequency at which one 
must apply pulses in order to maintain coherence. 

Figure [21 demonstrates the accuracy of the continuum 



GaAs Echo Decay Times 




Fock-Darwin radius [nm] 

FIG. 2: (Color online) Hahn echo total decay times, = 2tq , 
for GaAs quantum dots with various Fock-Darwin radii ro 
and thickness zq comparing the continuum approximation re- 
sults (solid lines) to 1 /e decay times obtained via the cluster 
expansion and previously presented in Ref. UlI This exam- 
ple includes dipolar but not indirect exchange interactions; it 
therefore slightly over-estimates the decay time. 

approximation compared with exact cluster calculations 
for Gaussian shaped quantum dots. It also shows that 
the continuum approximation is best for large dots and 
deviates as we consider smaller dots. In fact, it has been 
clearly reasoned 6 that the decay time must approach in- 
finity as the quantum dot size approaches zero extent, 
but the continuum approximation fails to capture this 
trend. 



D. When is the short-time limit appropriate? 

The short-time behavior will be exhibited on a time 
scale that is short relative to the time scale of the dy- 
namics of the relevant cluster contributions in . Be- 
cause a large number of these cluster contributions are 

added together in Re |(IFW) j [and then exponentiated 

to yield ve from Eq. (f2Tj)) ], it is possible for the decay 
of the echo, Ve, to occur on a time scale that is small 
relative to the dynamical time scales of any significantly 
contributing cluster. In particular, the decay will exhibit 
a short-time behavior when the clusters with the fastest 
dynamics dominate Re|(W^W)|. When there is a mix- 
ture of dynamical time scales playing a role, then the 
short-time behavior will be washed out by oscillations 
generated by HF-induced precessions of the nuclei. 

In considering the dynamical time scale of a cluster 
contribution, we really want to know the effect of this 
cluster on electron spin dephasing. This is determined 
by the difference in HF energies (with a reciprocal re- 
lationship between time and energy) for different spin 
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TABLE II: GaAs and/or InAs material fit parameters [see Eqs. (|46|l and (|47|) ] in microseconds with the magnetic field, B, 
in the [001] or [110] lattice directions. Considers only dipolar interactions except for values in parentheses that also include 
indirect exchange. 





B | 


[001] 




B || [110] 




Level 


M[100],[010] 


M[ooi] 


M[no] 


MfiToi 


M[ooi] 






GaAs with natural isotope abundances 




I = 


0.36 (0.29) 


0.33 (0.37) 


0.36 (0.28) 


0.28 (0.31) 


0.41 (0.30) 


I = 1 


0.25 (0.21) 


0.23 (0.26) 


0.26 (0.22) 


0.20 (0.20) 


0.29 (0.22) 


I = 2 


2.2 (1.8) 


2.0 (2.1) 


2.1 (1.7) 


1.7 (1.8) 


2.2 (1.8) 


I = 3 


4.1 (3.6) 


3.6 (3.7) 


3.8 (3.2) 


3.2 (3.3) 


3.9 (3.2) 






InAs with 100% 115 In 






I = 


0.29 


0.26 


0.30 


0.23 


0.34 


I = 1 


0.20 


0.19 


0.21 


0.16 


0.24 


I = 2 


1.6 


1.4 


1.5 


1.2 


1.6 


I = 3 


2.7 


2.4 


2.6 


2.2 


2.6 




As only (upper bound for GaAs/InAs mixtures) 




I = 


0.49-0.55 


0.45-0.50 


0.50-0.56 


0.38-0.43 


0.57-0.63 


I = 1 


0.35-0.39 


0.32-0.35 


0.36-0.40 


0.27-0.30 


0.40-0.45 


/ = 2 


3.3-3.8 


2.9-3.4 


3.2-3.6 


2.6-2.9 


3.2-3.7 


I = 3 


6.5-7.5 


5.6-6.5 


6.2-7.1 


5.1-6.0 


5.8-7.1 



polarization configurations of the nuclei in the cluster. 
In the extreme case that all of the HF energies are the 
same, there is no spectral diffusion induced (the dynami- 
cal time scale is infinite). The dynamical time scale is also 
determined by the interactions between the nuclei that 
can cause changes in the spin polarization configurations 
(turning these interactions off will also shut off spectral 
diffusion); however, we can estimate a lower bound time 
scale from just the inverse of differences in HF energies 
among the cluster. 

We first address the £l e — ► oo limit and later discuss, 
briefly, the short-time behavior of HF-mediated interac- 
tions. Our previous ^ 11 ' 12 results in the fl e — > oo limit 
show that the echo decay typically exhibits a short-time 
behavior in quantum dots with assumed Gaussian shaped 
wavefunctions but not in donors with exponential-like 
wavefunctions. This is understood in the following way. 
In the case of the donor, the fastest dynamics come from 
those few nuclei in the center that have large differences 
in HF energies. These are too few in number to domi- 
nate the decay; therefore, a mixture of time scales must 
play a role, slowing down but contributing more as we 
consider clusters further from the center, and the short- 
time behavior is washed out. For a quantum dot with a 
Gaussian shaped wavefunction, however, the fastest con- 
tributions occur where the wavefunction gradient is large 
in a ring around the dot at a radius of the characteristic 
size of the dot. There are many such clusters that are 
collectively capable of dominating the echo decay so that 
it will exhibit a short-time behavior. These arguments 
are illustrated in Figs. [3] and [4] 

There is a simple self-consistent check of the short- 
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FIG. 3: (Color online) A comparison of the position depen- 
dence of the electron probability density, absolute value of 
the gradient of this density [strictly speaking, || J^f^H and 
||J|P(r)|| for the top and bottom graphs respectively], and 
the resulting contribution to the short-time behavior of SD 
[integrand of Eq. (J46J) oc r (^P(r)) 2 and oc (£P(r)) 2 for the 
top and bottom graphs respectively] for a laterally confined 
quantum dot of the form of Eq. (|48p . The curves have ar- 
bitrary vertical scales. The short-time behavior exhibited by 
the SD decay of such quantum dots is related to the fact that 
the maximum contributions occur roughly where the gradient 
is maximum and differences in HF coupling among neighbor- 
ing nuclei are large. 
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FIG. 4: (Color online) A comparison of the position depen- 
dence of the electron probability density, absolute value of 
the gradient of this density, ||J^-P(r)|| oc P(r), and the re- 
sulting contribution to the short time behavior SD [integrand 
of Eq. (|46|l oc r 2 (-§pP{r)) ] for an electron with an exponen- 
tial probability density of the form P(f) = exp(— r/ro) where 
r = ||f|[. Donor bound electrons, such as in Si:P, have an 
exponential- type decay (though not quite as simple). The 
curves have arbitrary vertical scales. The failure of the time 
expansions in such systems is related to the fact that the 
maximum contributions do not occur where the gradient is 
maximum. 



time behavior [Eqs. (|45f and (|46]) ] in the Sl e — * oo limit 
by comparing t to the fastest dynamical timescale es- 
timated by the maximum gradient multiplied by the 
lattice spacing (as a typical distance scale between nu- 
clei). Thus, the short-time behavior is valid when 

To « 1/ <zmax|VA(f)| . For an electron probabil- 



ity density of the form of Eq. I|48p. max|V^4(r)| ~ 
max (A n /ro, A n /zo). With r$ 3> zq, the dynamical time 
scale has more to do with the position of the cluster in the 
z direction rather than the radial direction with a time 
scale estimate of t ~ Zo/aA n . With a 10 nm quantum 
dot thickness for zq and lattice constant of about 0.5 nm, 
this sets a time scale of about 10 fis. This estimation 
is slightly pessimistic because we have computed cluster 
expansion results for zq = 10 nm and ro = 100 nm to 
obtain To ~ 25 /is Hahn echo decay exhibiting the short- 
time behavior even though To > 10 /is. However, our 
10 ^s estimate is expected to be overly short because the 
dynamics is slowed by the weak 10 ms dipolar coupling 
(necessarily involved in any decoherence in the f2 e — * oo 
limit). 

This paper is primarily restricted to the £l e — > oo 
regime. For the moment, however, let us consider a fi- 
nite f2 e and discuss the question of short-time behavior 
validity with HF-mediated interactions. Because these 
are nonlocal interactions, differences in HF energies can 
be as large as the HF energies themselves. This sets a 
dynamical time scale on the order of 1/A n ~ 1 /xs. Also, 



the fact that HF-mediated interactions are weak does not 
help too much in this case because there are N > 10 5 nu- 
clei with which a given nucleus may interact; this yields 
a strong collective interaction of about NA nm ~ A n . In 
a manner of thinking, on a microsecond time scale, a 
given nucleus is likely to flip-flop with some other nu- 
cleus in the bath through the HF-mediated interaction. 
On the other hand, flip-flops between nuclei with large 
differences in HF energies will be suppressed due to HF 
energy conservation. For this reason, clusters with the 
fastest dynamics, those with large HF energy differences, 
will give weak contributions due to energy conservation, 
while the larger contributors with matching HF energies 
will have slow dynamics. This leads to a situation similar 
to that of the donor electron in the Q e — > oo limit, where 
the fastest contributors cannot dominate and, therefore, 
the short-time behavior is washed out. A proper treat- 
ment of HF-mediated interactions, therefore, would not 
use the time perturbation; instead, we should use the in- 
trabath perturbation, perturbatively treating both dipo- 
lar and HF-mediated interactions with respect to the HF 
interaction. 

The short-time behavior that we consider in this paper, 
with its convenient application in the continuum approx- 
imation, only applies in the £! e — > oo limit, where the 
HF-mediated interactions are negligible. How large does 
O e need to be for this limit to apply? We may pose this 
question differently to ask: How does fi e affect the time 
scale of short-time behavior convergence? If we can push 
this time scale sufficiently larger than To, the short-time 
behavior will emerge. The HF-mediated interaction rep- 
resents the lowest-order interaction with a magnetic-field 
dependence, A nm = A n A m /AQ, e . The time scale of the 
fast dynamics due to the HF-mediated interactions will 
be a combination of the difference in HF energies (~ A n ) 
and the collective HF-mediated energy (~ NA nm ) to 
give a time scale of ~ (A n NA nm ) 1 oc vOe- There- 
fore, according to this simple argument, it is necessary 
to quadruple fl e in order to double the time scale of 
the short-time behavior as long as HF-mediated interac- 
tions are dominant over local internuclear interactions. 
The lowest-order effects of HF-mediated interactions, 
however, are reversed by any DD refocusing technique. 
Therefore, this argument is only relevant for FID. 



SINGLET-TRIPLET DOUBLE QUANTUM 
DOTS 



Remarkable experiments have recently 3 -^ investigated 
the coherence properties of a single qubit in GaAs quan- 
tum dots. In the earlier of these experiments, 3 the qubit 
was not the spin of a single electron but, rather, a sub- 
space of two electron spins, each in separate quantum 
dots with a controllable exchange interaction between 
the two dots. The qubit states are represented by the 
two-electron spin states with zero total spin, ® |J.) 2 
and ® |t)2, where the 1 and 2 subscripts label the 
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dots (and contained electrons). An applied magnetic 
field protects each electron spin from depolarization; at 
the same time, the degeneracy of the zero-spin subspace 
is protected from uniform magnetic-field fluctuations. 45 
Electrostatic potentials are used to manipulate the elec- 
trons. State preparation and final readout are per- 
formed by biasing the two electrons, with an applied 
voltage, into the same dot so that the singlet state, 
(| T> 1 ® U>2 + ® Ith) /V2, has the lowest energy be- 
cause of the Pauli-exchange interaction^ 5 -^ Voltage con- 
trol is also used to turn on an exchange interaction by 
allowing the wavefunction of the two electrons on differ- 
ent dots to overlap; such control can be used to rotate 
the qubit i 45 ! 46 By using this control, one can apply ir 
pulses in order to perform a Hahn echo sequence or any 
other DD sequence (such as those discussed in Sec. Ill C| 
to prolong the coherence of the qubit. 

We can simply map this two-electron qubit into our 
single-spin qubit formalism. For convenience, we will 
define |0) = Jt)i ® ||) 2 and |1) = |4,> x ® |T>a as our 
two qubit basis states. Turning on the exchange inter- 
action will split the energies of the (|0) + |l|)/\/2 and 
(|0) — |1|)/a/2 superposition states and thereby rotate 
the qubit in a "transverse" direction as required for a 
DD sequence that combats dephasing. In order to ob- 
tain the free evolution Hamiltonian needed by our for- 
malism, we simply need to derive the qubit-bath Hamil- 
tonian, TCebS z , from the qubit-bath interactions in each 
of the two dots, H^Si z + S 2z , by taking its ma- 
trix elements in terms of our qubit basis states. With 
these definitions where we only have a dephasing cou- 
pling between the qubit and the bath, it is clear that 

(Q\H^S lz +H^S 2z \l) = (l\H^S lz +H^S 2z \0) = 0; 
we thus have only the following dephasing qubit-bath in- 
teraction: 

H eb = 2(0\H$S lz +H < $S 2z \0) 

= -2{l\n^S lz +H^S 2z \l) (53a) 
= H%-H$. (53b) 

For each dot i, the qubit-bath interaction is given by 
7Y W = - V A {l) I +- V A {l) I _J (5A) 

n£dot i n^m£dot i 

During the free evolution part of the pulse sequence, the 
two electrons must essentially have no overlap in their 
wavef unctions; therefore, will only be nonzero when 
n represents a nucleus in dot i. This is the justification 
for summing over only the relevant dot in Eq. (|54p . 

Assuming that the internuclear interactions occur only 
within the same bath (and that the bath is initially uncor- 
rected), then the problem fully decouples into spectral 
diffusion problems for dots l and 2 separately. With re- 
gards to in the cluster approximation [Eq. HO}], we 
simply need to sum the cluster contributions in the two 
dots separately. In a random unpolarized bath with two 



equivalent dots, the cluster contributions in each dot will 
be identical; then ve is simply the squared value of the 
echo for the problem of a single electron in just one of the 
dots. There should, thus, be no qualitative difference be- 
tween the spectral diffusion of a single-spin qubit and this 
double-spin qubit; a prediction of ve ~ exp [— (t/to) 4 ] 
for a single-spin qubit will carry over to the double-spin 
qubit. 

Although the reported Hahn echo decay time T 2 of 
Ref. is compatible with our theory (which disregards 
other decoherence mechanism) as a limiting case, it is 
clear that the experimental echo decay does not match 
the exp [— (t/to) 4 ] form. The experimentalists seem to 
be observing a decoherence mechanism that we are not 
treating. They report that the T 2 time increases with 
an increase in magnetic field^; therefore, they must not 
be operating in the high field limit regime. Our results 
may be viewed as yielding the best decoherence times 
achievable by increasing the applied magnetic field. 



A. Dynamic nuclear polarization and the 
"Zamboni effect" 

To minimize the effects of decoherence due to a bath of 
nuclear spins, one strategy is to polarize the nuclei. When 
they are polarized, they cannot flip-flop. This is particu- 
larly appealing in III-V semiconductors, where all of the 
isotopes have nonzero spin. Recent experiments have suc- 
cessfully achieved some degree of nuclear polarization in 
double quantum dot systems 4£ This is accomplished by 
biasing to a point where there is an anticrossing between 
the single state and the m s = +1 triplet state; the tran- 
sition between these states requires a nuclear spin flip to 
conserve angular momentum. By cycling through this 
anticrossing, they are able to produce polarizations of a 
few percent (producing effective nuclear fields of about 
20 mT in dots where full polarization would yield about 
5 T)*ffl 

Even with such modest polarization, there can be a 
significant impact on inhomogeneous broadening. It does 
this by effectively smoothing out the hyperfine field and, 
because of this smoothing, has been coined the Zam- 
boni effect by experimentalists^ Essentially, the process 
of dynamic nuclear polarization will most likely polarize 
those nuclei with the strongest coupling to the electron, 
those with the largest hyperfine coupling. These nuclei 
are also the most important in terms of inhomogeneous 
broadening (they give the largest contribution to the ef- 
fective magnetic field). Strong polarization is not nec- 
essarily involved in the suppression of inhomogeneous 
broadening. Homogenizing the system to remove the 
broadening need only change the polarization by roughly 
the same amount as the unbiased statistical broadening, 
which scales as 1/y/N; this is less than 1% for N = 10 5 . 
By removing the effects of inhomogeneous broadening in 
this way, it may be possible to view FID due to SD. 

This modest polarization will have a weak effect on 
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SD according to our theory. While T 2 (from inhomoge- 
neous broadening) is improved by this strategy, T 2 (from 
spectral diffusion) is not significantly altered. There are 
two reasons for this. First, the nuclei being polarized 
are not necessarily those nuclei responsible for SD. This 
is illustrated in Fig. [3] where the regions of electron oc- 
cupation probability do not correspond to the greatest 
SD contributors. Second, SD has a weak dependence 
on polarization because its contributors are clusters of 
two or more nuclei. Where we quantify polarization as 
p = P\ — P\ (the difference of the probability of being 
up versus down assuming spin 1/2 nuclei), the number 
of pairs that can flip-flop scales as (1 — p 2 )4£ When the 
spin is larger than 1/2, the dependence is even weaker 
(there is a larger fraction of states of two nuclei that can 
flip-flop). Therefore, one needs nearly 100% polarization 
in order to suppress SD T 2 decay. 

Therefore, we predict, with substantial confidence, 
that the coherence enhancement by the Zamboni effect 
could at best lead to a decoherence time of Tj(^> T 2 ), 
but never up to T 2 (> 7j), i.e., the Zamboni effect would 
never produce a coherence time longer than the spin echo 
coherence time. 



VI. DISCUSSION AND CONCLUSION 

The main result of this paper is the decoherence 
time formula of Eq. (|46|) with wavefunction dependence 
[Eq. (|52[) for the specific case of Gaussian-type simple 
harmonic oscillator confinement]. By using a table (Ta- 
ble|TT]for GaAs and InAs) of values for a time quantity we 

denote by /i, this formula will yield Tq^ for a given pulse 
sequence of I concatenations of the Hahn echo (1 = for 
free induction decay, and I = 1 for the Hahn spin echo 
itself); the initial behavior of the coherence or echo de- 
cay is exp(-[r/r^ ) ]P( / )) [Eqs. gT]) and ([51])]. where r is 
the time between pulses for the given sequence. By us- 
ing our definitions of decoherence times, Tj = rj^ 1 (free 

induction decay) and T 2 = (the traditional Hahn 

echo). The generalized concatenated echo decoherence 

times are T% = 2 1 Tq \ An interesting experimental test 
of this theory could be to compare the decoherence times 
of the same sample with an applied magnetic field in dif- 
ferent lattice directions and check for agreement with Ta- 
ble |TT] (for GaAs or InAs); in such a test, however, one 
must carefully account for any change in confinement as 
a result of changes in the applied magnetic field. For 
experimental systems that allow for the application of 
pulse sequences, with the ability to perform rapid ir ro- 
tations of the electron spin relative to dynamical time 
scales, a more straightforward test would be to compare 
different levels of concatenation and check for agreement 
with Table [Hand Eqs. dSJ) and gTJ). 

There are two important approximations that our de- 
coherence time formula assumes. First, we take the limit 
of a large applied magnetic field. This may not always 



be experimentally accessible, but, in any case, our re- 
sults represent the maximum coherence times achiev- 
able by applying a strong magnetic field. Second, we 
use a short-time approximation and discuss its valid- 
ity in Sec. IIVDI Failure of the short-time approxima- 
tion does not invalidate our general cluster expansion 
[Sec. IIII Al and Ref. [llj, however; it only means that our 
simple wavefunction dependent decoherence time formula 
[Eq. (|46|) ] is no longer accurate. 

Finally, in our considerations of the singlet-triplet (two 
electron) double quantum dot scenario, we show that it 
is equivalent to the single dot (one electron) case in terms 
of spectral diffusion assuming negligible exchange inter- 
action between pulses. (For a treatment that includes the 
exchange interaction, see Ref. H^.) We also predict that 
the Zamboni effect will have little impact upon spectral 
diffusion. 
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APPENDIX A: FACTORIZ ABILITY OF 
CLUSTER CONTRIBUTIONS 

We define W' c as the contribution of cluster C to W; 
that is, W' c gives the sum of all terms of W'- 1 ' that involve 
cluster C. Thus, 



E 

C^0 



W' e . 



(Al) 



Likewise, we will define W' A B as the sum of all terms of 
that involve clusters A and B so that 



WW = 



- E 

2 ^ 

Af\B=(D 



(A2) 



where the factor of 1/2 is necessary to compensate for 
the double counting of W' B A = W' A B and defining Ws 

to be the solution of W when only including nuclei in 
the set S (with all of the interactions between them). 
Similarly, we may define Wa,b to be the solution of W 
when only including nuclei in the sets A and B with 
interactions among A and among B but not between 
A and B. Because Wa,b is just a product of evo- 



lution operators of the form exp (- 



TL A + TL 
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Note 



exp (-il^r) exp (-iU%T^ , W AlB = W A ® W B . 

that W' c is the C cluster contribution to any Ws with 
S D C; this is simply due to the fact that any interac- 
tions of W that are not contained in Ws are irrelevant 
when considering terms that do not involve those interac- 
tions. For this reason, W' A is the A cluster contribution 
01 W a ,b — W A <g> Ws and W' B is its B cluster contribu- 
tion. Therefore, W' AB — W' A <£> W' B so that the double 
cluster contribution is simply the product of the individ- 
ual cluster contributions. 

This procedure may be applied to terms of any number 
of clusters so that 



{d } disjoint, <g> i 
Ci # 0: 



(A3) 



Assuming that the initial bath states are uncorrelated, 
then (W' A <g> W B ) = {W' A ) x (W B ) and 



(w) = e n^ 4 ), 



(A4) 



{Ci} disjoint, i 
Ci # 0, 



which essentially reproduces our cluster decomposition 
from Ref . [ll|. 

When only small clusters, relative to the size of the 
bath, give non-negligible contributions to W^ , then this 
factorizability allows us to make the following approxi- 
mation: 



(w [k] ) m ^(w [i] y 

kl 



(A5) 



assuming k <C N. The right hand side will provide all 
necessary products of cluster contributions without over- 



counting (the l/k\ factor compensates for permutation 
overcounting); however, it also includes products among 
overlapping clusters. In a large bath, sets of overlapping 
clusters are negligible compared to the number of sets of 
nonoverlapping clusters so that these extraneous terms 
are negligible. 



APPENDIX B: COMPUTING THE CONTINUUM 
APPROXIMATION TENSOR 



Computing MW of Eq. (43]) by calculating fjp(r m ) 
from Eq. (f3"4")l offers the advantage of simplifications 
due to the fact that operators which act on different 
sets of nuclei must commute. For example, to compute 

Inz , Hb , one need only consider the terms in Hj, that 

involve m. For us, however, it was more convenient to 
reuse a code that computes (W^) for any set of nuclei, C. 

By noting Eqs. (|21b|) and (|3"0|) , we may compute /„,„ by 
letting Ak oc 5 n ^ and summing together the lowest-order 
(in the time expansion) results of all cluster contribu- 
tions, (Wc), that include nucleus n (n e C). Similarly, 
if we let Ak oc S n ,k + &m,k for a given pair n ^ m, we 

+ fn% + fm,n- Subtracting off 
parts that may be computed by using 
Ak oc 5 n ,k and A k oc 6 m ,k, we can obtain + fm,n, 
which may then be used to compute M' 1 ' from Eqs. (|43|) 
and (|42|) . We may also use a statistical sampling of clus- 
ters to speed up the calculation of M>^, which is partic- 
ularly useful as one increases the number of concatena- 
tions, I. 



may compute f$ n + fm] m 
the fn'n and .in 



1 A. M. Tyryshkin and S. A. Lyon (private communication); 
A. M. Tyryshkin, S. A. Lyon, A.V. Astashkin, and A. 
M. Raitsimring, Phys. Rev. B 68, 193207 (2003); A. M. 
Tyryshkin, J. J. L. Morton, S. C. Benjamin, A. Ardavan, 
G. A. D. Briggs, J. W. Ager, and S. A. Lyon, J. Phys.: 
Condens. Matter 18, S783 (2006). 

2 E. Abe, K.M. Itoh, J. Isoya and S. Yamasaki, Phys. Rev. 
B 70, 033204 (2004); E. Abe, A. Fujimoto, J. Isoy a, S. 
Yamasaki and K.M. Itoh, arXiv:cond-mat/0512404 (un- 
published) . 

3 J. R. Petta (private communication); JR. Petta, A. C. 
Johnson, J. M. Taylor, E. A. Laird, A. Yacoby, M. D. 
Lukin, C. M. Marcus, M. P. Hanson, and A. C. Gossard, 
Science 309, 2180 (2005). 

4 F. H. L. Koppens, J. A. Folk, J. M. Elzerman, R. Han- 
son, L. H. Willems van Beveren, I. T. Vink, H. P. Tranitz, 
W. Wegscheider, L. P. Kouwenhoven, and L. M. K. Van- 
dersypen, Science 309, 1346 (2005); R. Hanson, L. H. 
Willems van Beveren, I. T. Vink, J. M. Elzerman, W. J. M. 
Naber, F. H. L. Koppens, L. P. Kouwenhoven, and L. M. 
K. Vandersypen, Phys. Rev. Lett. 94, 196802 (2005); F. 



H. L. Koppens, K. C. Nowack, and L. M. K. Vandersypen, 
larXiv:0711.0479l (unpublished). 

5 A. Greilich, D. R. Yakovlev, A. Shabaev, Al. L. Efros, I. A. 
Yugova, R. Oulton, V. Stavarache, D. Reuter, A. Wieck, 
and M. Bayer, Science 313, 341 (2006). 

6 R. de Sousa and S. Das Sarma, Phys. Rev. B 68, 115322 
(2003). 

7 D. Klauser, W. A. Coish, and D. Loss, Adv. Solid State 
Phys. 46, 17 (2007). W. A. Coish and D. Loss, Phys. Rev. 
B 70, 195340 (2004). 

8 W. M. Witzel, R. de Sousa, and S. Das Sarma, Phys. Rev. 
B 72, 161306(R) (2005). 

9 Wang Yao, Ren-Bao Liu, and L. J. Sham, Phys. Rev. B 
74, 195301 (2006). 

C. Deng and X. Hu, Phy s. Rev. B 73 241303(R) (2006); 
arXiv:cond-mat /0608544 (unpublished) . 

1 W. M. Witzel and S. Das Sarma, Phys. Rev. B 74, 035322 
(2006). 

2 W. M. Witzel and S. Das Sarma, Phys. Rev. Lett. 98, 
077601 (2007). 

3 Wang Yao, Ren-Bao Liu, and L. J. Sham, Phys. Rev. Lett. 



17 



14 



98, 077602 (2007). 



Ren-Bao Liu, Wang Yao, and L. J. Sham, New J. Phys. 9, 
226 (2007). 

15 W. M. Witzel and S. Das Sarma, Phys. Rev. B 76, 
241303(R) (2007). 

16 B. Lee, W. M. Witzel, and S. Das Sarma. larXTv:0710.1416l 
(accepted to Phys. Rev. Lett.). 

17 Gotz S. Uhrig, Phys. Rev. Lett. 98, 100504 (2007). 

18 B. Herzog and E. L. Hahn, Phys. Rev. 103, 148 (1956); A. 
M. Portis, ibid. 104, 584 (1956). 

19 G. Feher and E. A. Gere, Phys. Rev. 114, 1245 (1959). 

20 W. B. Mims and K. Nassau, Bull. Am. Phys. Soc. 5, 419 
(1960); W. B. Mims, K. Nassau, and J. D. McGee, Phys. 
Rev. 123, 2059 (1961). 

21 J. R. Klauder and P. W. Anderson, Phys. Rev. 125, 912 
(1962). 

22 G. M. Zhidomirov and K. M. Salikhov, Sov. Phys. JETP 
29, 1037 (1969). 

23 M. Chiba and A. Hirai, J. Phys. Soc. Jpn. 33, 730 (1972). 

24 S. K. Saikin, Wang Yao, and L. J. Sham, Phys. Rev. B 75, 
125314 (2007). 

25 R. de Sousa and S. Das Sarma, Phys. Rev. B 67, 033301 
(2003). 

26 R. F. Shulman, J. M. Mays, and D. W. McCall, Phys. Rev. 
100, 692 (1955). 

27 R. F. Shulman, B. J. Wyluda, and H. J. Hrostowski, Phys. 
Rev. 109, 808 (1958). 

28 R. K. Sundfors, Phys. Rev. 185, 458 (1969). 

29 W. M. Witzel, Xuedong Hu, S. Das Sarma, Phys. Rev. B 
76, 035212 (2007). 

30 R. Winkler, Spin- Orbit Coupling Effects in Two- 
Dimensional Electron and Hole Systems (Springer-Verlag, 
Berlin, 2003), Appendix B. 

31 C. P. Slichter, Principles of Magnetic Resonance, 3rd ed. 
(Springer-Verlag, Berlin, 1990). 

32 N. Bloembergen and T. J. Rowland, Phys. Rev. 97, 1679 
(1955). 

33 N. Shenvi, R. de Sousa, and K. B. Whaley, Phys. Rev. B 
71, 224411 (2005). 



34 A. Abragam, The Principles of Nuclear Magnetism (Ox- 
ford University Press, London, 1961), Chap. IV, Eq. (63). 

35 P. W. Anderson, Phys. Rev. 99, 623 (1955). 

36 K. Khodjasteh and D. A. Lidar, Phys. Rev. A. 75, 062310 
(2007). 

37 W. Magnus, Commun. Pure Appl. Math. 7, 649 (1954). 

38 L. Viola and E. Knill, Phys. Rev. Lett. 94, 060502 (2005). 

39 K. Khodjasteh and D.A. Lidar, Phys. Rev. Lett. 95, 
180501 (2005). 

40 H. Y. Carr and E. M. Purcell, Phys. Rev. 94, 630 (1954); S. 
Meiboom and D. Gill, Rev. Sci. Instrum. 29, 6881 (1958). 

41 C. Kittel, Introduction to Solid State Physics, 7th ed. (Wi- 
ley, New York, 1996). 

42 CRC Handbook of Chemistry and Physics, 70th ed., E-82 
(CRC, Boca Raton, FL, 1989). 

43 Encyclopedia of Nuclear Magnetic Resonance, D. M. 
Granty and R. K. Harris, (eds.), vol. 5, John Wiley & 
Sons, Chichester, UK, 1996. 

44 D. Paget, G. Lampel, B. Sapoval, and V. I. Safarov, Phys. 
Rev. B 15, 5780 (1977). 

45 J. M. Taylor, H.-A. Engel, W. Dur, A. Yacoby, C. M. Mar- 
cus, P. Zoller, and M. D. Lukin, Nat. Phys. 1, 177 (2005). 

46 J. R. Petta, A. C. Johnson, J. M. Taylor, A. Yacoby, M. D. 
Lukin, C. M. Marcus, M. P. Hanson, and A. C. Gossard, 
Physica E (Amsterdam) 34, 42 (2006); E. A. Laird, J. R. 
Petta, A. C. Johnson, C. M. Marcus, A. Yacoby, M. P. 
Hanson, and A. C. Gossard, Phys. Rev. Lett. 97, 056801 
(2006); J. M. Taylor, J. R. Petta, A. C. Johnson, A. Ya- 
coby, C. M. Marcus, and M. D. Lukin, Phys. Rev. B 76, 
035315 (2007). 

47 J. R. Petta, J. M. Taylor, A. C. Johnson, A. Yacoby, M. 
D. Lukin, C. M. Marcus, M. P. Hanson, and A. C. Gos- 
sard, Phys. Rev. Lett. 100, 067601 (2008); C. M. Marcus 
(private communication) . 

48 S. Das Sarma, Rogerio de Sousa, Xuedong Hu, and Belita 
Koiller, Solid State Commun. 133, 737 (2005). 

49 W. Yang and R. B. Liu, Phys. Rev. B 77, 085302 (2008). 



